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Recent work has highlighted remarkable effects of classical thermal fluctuations in the dipo- 
lar spin ice compounds, such as "artificial magnetostatics", manifesting as Coulombic power- 
law spin correlations and particles behaving as diffusive "magnetic monopoles". In this pa- 
per, we address quantum spin ice, giving a unifying framework for the study of magnetism 
of a large class of magnetic compounds with the pyrochlore structure, and in particular dis- 
cuss Yb 2 Ti 2 7 and extract its full set of Hamiltonian parameters from high field inelastic 
neutron scattering experiments. We show that fluctuations in Yb 2 Ti 2 7 are strong, and that 
the Hamiltonian may support a Coulombic "Quantum Spin Liquid" ground state in low field 
and host an unusual quantum critical point at larger fields. This appears consistent with 
puzzling features in prior experiments on Yb 2 Ti 2 7 . Thus Yb 2 Ti 2 7 is the first quantum 
spin liquid candidate in which the Hamiltonian is quantitatively known. 

Rare earth pyrochlores display a diverse set of fascinating physical phenomena. 1 One of the 
most interesting aspects of these materials from the point of view of fundamental physics is the 
strong frustration experienced by coupled magnetic moments on this lattice. The best explored ma- 
terials exhibiting this frustration are the "spin-ice" compounds, Ho 2 Ti 2 7 , Dy 2 Ti 2 7 , in which the 
moments can be regarded as classical spins with a strong easy-axis (Ising) anisotropy. 2 3 The frus- 
tration of these moments results in a remarkable classical spin liquid regime displaying Coulombic 
correlations and emergent "magnetic monopole" excitations that have now been studied exten- 
sively in theory and experimentP^ 

Strong quantum effects are absent in the spin ice compounds, but can be significant in other 
rare earth pyrochlores. In particular, in many materials the low energy spin dynamics may be re- 



duced to that of an effective spin S = 1/2 moment, with the strongest possible quantum effects 
expected. In this case symmetry considerations reduce the exchange constant phase space of the 
nearest neighbour exchange Hamiltonian to a maximum of three dimensionless parameters. 7 The 
compounds Yb 2 Ti 2 7 , Er 2 Ti 2 7 , Pr 2 Sn 2 7 ! (and possibly Tb 2 Ti 2 07®) are of this type, and it has 
recently been argued that the spins in Yb 2 Ti 2 7 and Er 2 Ti 2 7 are controlled by exchange coupling 
rather than by the long-range dipolar interactions which dominate in spin iceJ2U2l This makes these 
materials beautiful examples of strongly quantum magnets on the highly frustrated pyrochlore lat- 
tice. They are also nearly ideal subjects for detailed experimental investigation, existing as they do 
in large high purity single crystals, and with large magnetic moments amenable to neutron scatter- 
ing studies. Yb 2 Ti 2 7 is particularly appealing because its lowest Kramers doublet is extremely 
well separated from the first excited one, 11 and a very large single crystal neutron scattering data 
set is available, allowing us to measure the full Hamiltonian quantitatively, as we will show. Al- 
though we specialize to the latter material in the present article the theoretical considerations and 
parameter determination method described here will very generally apply to all pyrochlore materi- 
als where exchange interactions dominate and whose dynamics can be described by that of a single 
doublet. 

Theoretical studies have pointed to the likelihood of unusual ground states of quantum antifer- 
romagnets on the pyrochlore lattice.E3i2 Most exciting is the possibility of a quantum spin liquid 
(QSL) state, which avoids magnetic ordering and freezing even at absolute zero temperature, and 
whose elementary excitations carry fractional quantum numbers and are decidedly different from 
spin wavespl Intriguingly, neutron scattering measurements have reported a lack of magnetic or- 



dering and the absence of spin waves in Yb 2 Ti 2 7 at low fields LL5H3 In a recent study, sharp spin 
waves emerged when a magnetic field of 0.5T or larger was applied, suggesting that the system 
transitioned into a conventional state. 16 The possible identification of the low field state with a 
quantum spin liquid is tantalizing, but progress certainly requires a more detailed understanding of 
the spin Hamiltonian. 

In this article, we present a detailed experimental and theoretical investigation of the excitation 
spectrum in the high field state throughout the Brillouin zone. We show that the spectrum is 
extremely well fit by spin wave theory, and through this fit we unambiguously extract all the 
microscopic exchange parameters in the spin Hamiltonian (see below). Interestingly, we find that 
the largest exchange interaction is of the same "frustrated ferromagnetic" Ising type as in spin 
ice, despite the fact that the g-tensor tends to orient magnetic moments primarily normal to the 
Ising axes. Moreover, spin-flip terms which induce quantum dynamics are comparable to the Ising 
exchange, which confirms the picture of Yb 2 Ti 2 07 as a strongly quantum magnet, in qualitative 
agreement with recent studies. ^ 17 ! 18 ! Strikingly, we find that the predictions of mean field theory 
(MFT) using these parameters disagree drastically with experiment in zero field, indicating that 
fluctuations strongly reduce or destroy the classically expected spin order. Taken together, these 
observations make Yb 2 Ti 2 7 a promising candidate for observation of QSL physics. The precise 
determination of the microscopic spin interactions sets the stage for a quantitative understanding 
and test of this proposal. 

Time- of- flight neutron scattering measurements were performed on a 7g single crystal of Yb 2 Ti 2 7 , 



grown via the optical floating zone method. Details of the crystal growth were given elsewhere ! 16 ! 19 ! 
The neutron scattering data was collected at the Disk Chopper Spectrometer at the NIST Center for 
Neutron Research, using 5 A incident neutrons. This configuration allowed an energy resolution 
of 0.09meV. The sample environment consisted of a 11. 5T vertical field magnet combined with a 
dilution insert that achieved a base temperature of 30mK. The scattering plane was HHL, with the 
field applied along the [110] vertical direction. The sample was rotated 147 degrees in 1.5 degree 
steps about the vertical, allowing a three dimensional data set to be acquired, i.e. two components 
of the wavevector Q within the scattering plane, and energy transfer. The spin excitation spectra 
along several high symmetry directions with the scattering plane were thereby obtained. 

At 30mK, the inelastic spectrum changes qualitatively at H = 0.5T; above this field strength, 
resolution-limited spin wave excitations which go soft with quadratic dispersion at nuclear-allowed 
positions develop, indicating a transition to a field-polarized ferromagnetic stated The spin wave 
excitations indicate that the symmetry of the underlying lattice is preserved, as is evident from 
gaps in the spectrum at the nuclear zone boundaries. In Figure 1 we show the spin wave disper- 
sions along several directions in the HHL plane for H = 2T and H = 5T These high symmetry 
directions are shown relative to the Brillouin zone structure within the HHL plane in Fig. 2. 

We compare the experimental data to spin-wave theory. We assume nearest- neighbour ex- 
change coupling only, as appropriate to the strongly localized /-electron states, and neglect dipolar 
interactions. The Hamiltonian, written in global spin coordinates, is then 

H = \ E jpsrS] - n B H» £ gTS?, (1) 

ij i 



where J?? = J^ is the matrix of exchange couplings between sites i and j, g^ v is the g-tensor 
for the spin at site i, and we take % — 1. Symmetry allows four independent exchange constants,^ 
Ji, • • • , J 4 . To specify them, we give the exchange matrix on one pair of nearest-neighbour sites, 
located at positions r = |(1, 1, 1) and r x = |(1, —1, —1) on a tetrahedron centered at the origin 

(a is the conventional cubic lattice spacing for the fee Bravais lattice): 

(J2 Ja J a \ 
-■h -h h\. (2) 

-h h h I 

The other exchange matrices can be obtained from this one by cubic rotations given in the Supple- 
mentary Information. The g-tensor contains two components: g z parallel to and g xy perpendicular 
to the local C 3 rotation axis through the Yb site. 

Spin wave theory, carried out as described in the Supplementary Information, is fit to the H = 
5 T, T = 30 mK measurements; the fitting procedure focuses on the dispersion relation alone, and 
the overall intensity of the calculated spin waves is scaled to agree with the experiment at a single 
wavevector and energy point. The resulting inelastic structure factor 5(Q, u) (see Supplementary 
Information) is compared to both the 5T and 2T data in Figure 1 . The best fit is achieved with the 
following exchange parameters, in meV: 

J- L = -0.09 ± 0.03, J 2 = -0.22 ± 0.03, J 3 = -0.29 ± 0.02, J 4 = 0.01 ± 0.02. (3) 

Here we quote rough uncertainties obtained by the visual comparison of the theoretical and exper- 
imental intensities. The fit is performed by taking the ratio of components of the g-tensor to be 



9xy/9z = 2.4, i.e. the ratio obtained by Ref. 11 The fit then produces g z = 1.80, in nearly perfect 
agreement with these studies (using the g-factor ratio of Cao et al. instead, 17 i.e. g xy /g z = 1.8, does 



not reproduce the data as precisely). Using these results, a high temperature expansion gives (see 
Supplementary Information) a theoretical Curie- Weiss temperature ®cw = 312mK, which is com- 
parable to but smaller than the experimentally determined values, Qcw = 400mKpland 750mKp^ 
The deviations may be explained by the sensitivity of the theoretical value to small changes in 
the g-factors and exchange parameters, and to the dependence of the experimental value on the 
fitting range. 17 Furthermore, and most importantly, our extracted exchange parameters correctly 
reproduce relative intensities as well as the shape of the spin wave dispersion for each of the five 
directions. Agreement is excellent for H = 2T, showing that these parameters produce a robust 
description of the field-induced ferromagnetic state. We note, however, that there is a significant 
quantitative disagreement with the exchange parameters quoted in Refs.^, 10 (see Supplementary 
Information). 

Implications: The excellent agreement with spin-wave theory for fields H > 2T clearly 
indicates that the high field state is accurately modeled semi-classically, and is smoothly connected 
to the fully polarized limit. Theoretically, the ground state in this regime breaks no symmetries, and 
supports a ferromagnetic polarization along the axis of the applied field (for the (1 10) field used in 
the experiment). However, the semiclassical analysis clearly and dramatically fails at small fields, 
where the measurements show no signs of spontaneous long range orderJ-SI The classical zero 
field ground state for our Hamiltonian parameters has a large spontaneous polarization along the 
(100) axis. Extension of this analysis to a T > mean- field theory wrongly predicts a continuous 
magnetic ordering transition at a temperature of T C MF = 3.2K (see Supplementary Information). 
The experimental indications of a zero-field transition to long range order are mixedpiE^but early 



specific heat measurements^ found an anomaly at T c = 214mK, and Mossbauer spectroscopy^ 
suggested a transition at 240mK. This temperature is approximately 14 times lower than T C MF . 
If there is magnetic ordering at all, it appears to be substantially suppressed, indicating strong 
fluctuations - classical, quantum, or both. 

The presence of strong fluctuations makes a QSL ground state plausible in low field. We now 
use the Hamiltonian parameters to suggest the nature of this state. To do so, we rewrite the zero 
field Hamiltonian in terms of spins quantized along the local C3 axis for each site, similarly to Ref. 

h = ^{j 2 A^-J ± (s+S7 + srs+) + j ±± [ 7j A + s; + 7^S7" 

+J Z± [SftCySt + QS") +i<*j]}, (4) 

where here Sf are local spin coordinates, J zz = — |(2 J x — J2 + 2( J3 + 2 J 4 )), J± = \{2J\ — J 2 — 
J 3 - 2J 4 ), J ±± = \{Ji + J 2 - 2J 3 + 2J 4 ) and J z ± = 7^{Ji + J2 + J3 - -h), and the matrices 
7jj, Qj consist of unimodular complex numbers (given in the Supplementary Information). From 
the fits in Eq. Q, we find, in meV, 

J zz = 0.17 ± 0.04, J± = 0.05 ± 0.01, J ±± = 0.05 ± 0.01, J z ± = -0.14 ± 0.01, (5) 

where the uncertainties have been estimated by treating those in Eq. ^ as Gaussian random vari- 
ables. Note that the strongest interaction is J zz > 0, which precisely coincides with the nearest- 
neighbour spin-ice model. The model with J± and J zz only has been studied theoretically ! 13 l 24 l 
It does indeed support a QSL ground state, for sufficiently small J±/J zz . For larger J±/J zz , the 
ground state is instead a magnet with (Sf) ^ 0.^ While the actual value of J±/J zz ~ 0.3 



would place this model in the magnetic stated the J z ± interaction in particular is non-negligible 
in Yb2Ti 2 07, and preliminary theoretical work suggests that it tends to stabilize the QSL state. 
Indeed, in perturbation theory, the leading effect of the J z ± coupling is to induce in the effective 



Hamiltonian a term close to the Rokhsar-Kivelson interaction of Ref. 13 (see Supplementary In- 
formation), which was shown to stabilize the QSL.LL^l Although perturbation theory is strictly 
speaking only valid for J z ±/J zz <C 1, the conclusions of this analysis are likely to extend to a 
larger range of values. A non-perturbative study of the full Hamiltonian in Eq. ([4]) is beyond the 
scope of this paper, but will be reported in a future publication. Given the uncertainty in the phase 
boundary for the QSL state, we cannot disregard the possibility that an ideal sample of Yb 2 Ti 2 07 
would be magnetically ordered, but that such order is here disrupted by crystal defects. This pos- 
sibility should be pursued further in the future taking into account the new understanding of the 
Hamiltonian. We proceed now to discuss the implications of the alternative possibility of a zero 
field QSL state. 



Many of the key properties of the QSL state of the J zz — J± model were established in Ref. [13 
Conceptually, it is a quantum analog of the classical regime of Coulombic correlations observed in 
spin iceP Specifically, where spin ice realizes an analog of magneto statics, the QSL state of Eq. Q 
embodies a complete fictitious quantum electrodynamics. In this phase, the magnetic monopoles 
of spin ice become full-fledged coherent excitations of the system. In addition, the QSL supports 
dual electric monopoles and a dynamical emergent photon mode at low energy. The complex 
and largely diffuse character of the scattering in zero fielcP^ may well be a consequence of the 
combination of these diverse excitations. Indeed, where a neutron can create just one spin wave 



(S = 1) excitation, the S = 1/2 magnetic monopoles are excited in pairs with no individual 
momentum conservation constraints. A careful comparison of theoretical modeling and focused 
experiments in low field is clearly needed. 

A key consequence of the QSL scenario is the presence of a quantum confinement phase tran- 
sition in applied field (see Fig. 3). Such a quantum phase transition is required to remove the 
"fractional" excitations of the QSL phase (electric and magnetic monopoles) from the spectrum in 
the semiclassical high field phase. Theoretically, such quantum critical points have been studied 
in related model Hamiltonians, and occur by a mechanism analogous to the Higgs transition in the 
standard modeP^ or by magnetic monopole condensation. 27 The gapless excitations observed by 



neutrons at H « 0.5T in Ref. 16 indeed indicate a quantum critical point at this field. Further 



theoretical work is required for detailed comparison to experiment. 

Yb 2 Ti 2 7 enjoys several major advantages over other materials considered candidates for QSL 
states up to now!-^ it is the only case in which the Hamiltonian parameters are precisely known, and 
for which large single-crystal samples highly suitable for detailed neutron scattering measurements 
are available. Moreover, similar methods may be applicable to other rare earth pyrochlores in 
which substantial quantum effects are present, such as Er 2 Ti 2 07 and Tb 2 Ti 2 7 .GE Especially in 
this broader context, the prospects for detailed observation of the long-sought QSL physics are 
bright. The basic framework established here will allow coordinated theoretical and experimental 
studies to confront the problem. 
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Figure 1: The measured S(Q, uS) at T = 30mK, sliced along various directions in the HHL 
plane, for both H = 5T (first row) and H = 2T (third row). The second and fourth row show the 
calculated spectrum for these two field strengths, based on an anisotropic exchange model with five 
free parameters (see text) that were extracted by fitting to the 5T data set. For a realistic comparison 
to the data, the calculated S(Q,ui) is convoluted with a gaussian of full-width 0.09meV. Both 
the 2T and 5T data sets, comprised of spin wave dispersions along five different directions, are 
described extremely well by the same parameters. 



22L 



HHO 



HH1 




[110] direction 

Figure 2: Representations of the HHL scattering plane, showing the FCC Brillouin zone bound- 
aries and the corresponding zone centers (labelled in terms of the conventional simple-cubic unit 
cell). Blue lines indicate the directions of the five different cuts shown in Figure 1. 
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Figure 3: Schematic phase diagram in the temperature (T) - magnetic field (H) plane, for a 
material in the quantum spin liquid (QSL) phase of Eq. (1) at T = H = 0. At low field and tem- 
perature, the QSL state supports exotic excitations: magnetic (red sphere) and "electric" (yellow 
sphere) monopoles, and an emergent photon (wavy line). The field H c marks a quantum critical 
point: the confinement phase transition. For H > H c , the ground state is a simple field-polarized 
ferromagnet (FM), and the elementary excitations are conventional magnetic dipoles. The grada- 
tions in gray scale indicate crossovers to the quantum critical region between them, governed by 
the T = confinement phase transition. 
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Supplementary Information 



1 Cubic Rotations and Local Bases 

As described in the main text, we use the usual coordinate system for the pyrochlore lattice, with 
sites located on tetrahedra whose centers form a FCC lattice. We take one to be centered at the 
origin with its four corners at r = |(1, 1, 1), r x = |(1, —1, —1), r 2 = |(— 1, 1, —1) and r 3 = 
|(— 1,— 1,1). The exchange matrices J^ between sites of types i and j are obtained by applying 
the following cubic rotations IZij to J i: 

• 7^02 is a ^ rotation about the [111] axis, 

• TZ 03 is a ^ rotation about the [111] axis, 

• 1Z 2 i is a ^ rotation about the [111] axis, 

• IZ31 is a ^ rotation about the [111] axis, 

• IZ23 is a rotation made of a ^ rotation about the [111] axis followed by a ^ rotation about 
the [111] axis. 



Note Jji = Jfj. 
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We use the following local (a», b i: e,) bases 



eo = 73(1,1,1) 
fii = ^(1,-1,-1) 

e 2 = ^(-1,1,-1) 

e 3 = ^(-1,-1,1), 



y/E^ "• " 



-2,1,1' 



a 
ai 
a 2 = ^(2,l,-l; 

I ^3 = ^(2,-1,1) 



Ta^-'- 1 



V6' 



(6) 



bj = &i x aj, and the 4x4 complex unimodular matrices 



( o 



-1 e*3 e l 3 \ 



-1 





— 7 — 

e s 


7 — 

e s 


7- 
e 3 


e 3 





-1 


Ve"^ 


i- 

e s 


-1 


/ 



7 = -C, 



(7) 



for which our exchange Hamiltonian takes the form of Eq. (4). 



2 Curie- Weiss Temperature 



A high-temperature expansion of the Hamiltonian of Eq. (1) yields the 0(1/T 2 ) term in the uni- 
form susceptibility from which we extract the Curie- Weiss temperature, 



e 



cw 



Qk B {2gl y + gl) 



[2g 2 x M + J 2 - 5J 3 + 2J 4 ) 



(8) 



+Sg xy g z (J 1 + J2 + J3- J a) + g 2 z {2Ji -J 2 + 2J 3 + 4J 4 )], 



where k B is the Boltzmann constant. Using the formulation of Eq. (4) Qcw takes the simpler form 



e 



cw 



2k B {2gl y + gl) 



[g 2 z J zz -4g 2 xy (J± + 2J±±)-8V2g xy g z J z ±]. 



(9) 



3 Spin Wave Theory 



As usual, we expand the Hamiltonian about one of the classical states using Holstein-Primakoff 
bosons in the spirit of large s, and keep only terms of and up to order s, which shall then be set to 
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1/2. We define the transverse Holstein-Primakoff bosons x a = x a , y a = y\, conjugate with one 
another on site a of the pyrochlore lattice, satisfying 

[Xa,Va]=i, n «= X -^Y^~\ < 10 ) 

such that 

S a -u a = s-n a , S - v a = y/sx a , S a -w a = ^fsy a . (11) 

Here (v a , w a ,u a ) is an orthonormal basis, chosen so that u a gives the direction of the spin in 
the classical ground state at site a (which we find numerically). We find that for all fields, the 
ground state does not enlarge the unit cell, so that there are only four distinct such bases, which 
we denote by a = 0, ..,3. One may choose, for example, v a = u a x (1, 1, l)/||u a x (1, 1, 1)|| and 

W a = U a X v a . 

Since the classical ground state does not enlarge the unit cell, we can readily proceed to Fourier 
space in the four site basis. Keeping only terms of order s, we arrive at the spin-wave (quadratic) 
Hamiltonian, 



H,={xi k r 4 j^f z)[y:)> (i2) 

where ( X T Y T ) = ( x .. x 3 y .. y 3 ) andDf = D ab cos(k-(r a -r b )), D = A,B,C, 
with 



Kb = sv a ■ 3 ab ■ v b + ^H • g a • u a 5 ab - su a ■ Y^ J ab • u b 

1 b 

Bab = SW a ■ 3 ab ■ W b + ^H • g a • U a 5 ab - SU a -^2 J a6 ' U 6 

1 b 

Cab = SV a • J ab ■ W b , (13) 
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where H is the magnetic field and 3 ab and g a are 3x3 matrices with matrix elements Jf£ and 
g% v , respectively. To find the modes, we resort to the path-integral formulation. The action at 
temperature T = l/(k B /3) is 

s = ^T,i: z \-* n [Gk + (^n)r] z^ n , (14) 

Z P n k 



where u n = ^p is the Matsubara frequency, where we have defined Z T = ( X T Y T J , 

^ = 2(4? %) and r=f, 1 ° ~i U ) (15) 



(1 4 is the four-by-four identity matrix). As usual, the real frequency dispersion relations w(k) are 
found by solving the zero eigenvalue equations of the matrix G k + coT . Here, these are equivalently 
the (both zero and non-zero) eigenvalues of — TG^. 

We calculate the inelastic structure factor (to which the intensity J(k, ui) of the scattering is 
proportional) obtained from the moment-moment correlation function,^ 

S(k,u) = X) [<W - (k) M (k),] E«(" k - -wK(k,w)>, (16) 

where k is the unit vector associated with k, m£(k, o>) = HBH K ,9a K ^a^ UJ ) i s me moment at 
momentum k and real frequency ui on the a = 0, .., 3 sublattice in direction /i = x,y,z, and S£ 
is, as usual, the Kth coordinate of the effective spin-1/2 spin at site a. J^^k) = S^ — (k) M (k) v 
selects the component of the spin-spin correlations perpendicular to the scattering vectorpS and 
Gabfe, of) = (mj(- k, — co)ml(k, to)) is the moment-moment correlation function that originates 
from the interaction between the neutrons' moment and the spins' moments, which we find to be 

QZ = -^«E^( W - e a ) I (bPl,aW R , a } b ) , (17) 
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where 5 is the Dirac distribution, r]% = J2 K =x, y ,z 9% K ( v 1 w a\ e a is the ath eigenvalue of — TG, 
and ipR t a is its corresponding "right" eigenvector, i.e. such that — TG ■ ipR ;a = e« ^R,a- Also note 
that the momentum and frequency dependences are implied everywhere to be k, u. 

Now we estimate the amplitude of quantum fluctuations using our spin wave theory, by eval- 
uating the quantum moment reduction in zero field and at zero temperature. In the Holstein- 
Primakoff boson language the reduction of the spin expectation value from the classical value 
of 1/2, averaged over the four site basis, is r = | Sa=o(0| n a|0)> where (0|n a |0) indicates a 
ground state quantum expectation value. Evaluating this using the path integral method, we obtain 
r=|/ k)U; Tr([G k + a;r]- 1 )-|,i.e. 

r= o[L e WT^ 9 ^0.05 (18) 

(0 is the usual Heaviside distribution), which is a 10% reduction compared with the classical value 

of 1/2. 

4 Further Details of Experimental Data and Fits 

The inelastic neutron scattering data (rows 1 and 3 in Figure 1) contains several features worth 
commenting on further. First, the darkest blue areas do not contain any data, either for kinematic 
reasons near Q = (0, 0, 0), or because of the finite angular extend of the scan. Second, at E = 
one observes intensity which is due to coherent and incoherent elastic scattering from the sam- 
ple, and hence is more intense than the inelastic features by up to several orders of magnitude, 
thus appearing red (off scale). Third, near the (0, 0, 0) position there is higher background lead- 
ing to unphysical intensities due to contamination from the un-scattered incident beam. This is 
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observable in the HHH and HHO data sets near the zero position. 

The data sets were collected by counting for 8 minutes per angular rotation of the sample. The 
total time per magnetic field setting was about 12 hours. 

Supplementary Figure 1 shows the dispersions obtained from the fitting procedure over-plotted 
on the data. These fits were accomplished by digitizing the shape of the dispersion from the 
experimental data, and performing a least squares minimization routine to match it. Intensities 
were calculated based on the spin wave theory using the extracted exchange parameters, and were 
not fit to the data. The 5D parameter space (four exchange plus one g-tensor parameter) was 
explored using a uniform search technique with the resulting excellent description obtained (Figure 
1). 






-0.5 0.5 1 1.5 -2 2-2 2 12 1 2 

HHH(r.l.u.) 11L(r.l.u.) 22L(r.l.u.) HHO (r.l.u.) HH1 (r.l.u.) 

Supplementary Figure 1: Dispersions (white curves) obtained from the fitting procedure over- 
plotted on the data at H = 5 T 



5 Features of the Spin Wave Spectrum 



The spectra of Figure 1 are very rich. Because of the many parameters and of the presence of 
classical phase transitions (there are several ground states to the classical model), it is very difficult 
to track which features are due to what terms of Eq. (4). One remarkable feature which can be 
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easily identified throughout most of phase space (including the region around which Eq. (3) places 
Yb 2 Ti 2 07) is a quasi-flat band. Specifically, one spin wave mode is completely dispersionless in 
the plane of reciprocal space with k x = k y , i.e. normal to the magnetic field direction and passing 
through the origin in reciprocal space. All scattering measurements on Yb 2 Ti 2 7 have been taken 
in this plane, so this feature is quite significant in the experiments. In the region of phase space 
around Yb 2 Ti 2 7 and for H = 5 T, we find its energy to be, numerically: 

£ 2 dfiat « 0.74 + 0.51 J zz - 1.18 J± - 3.11 J±± - 5.81 J z ± meV, (19) 

where the J a 's must be input in meV, and which for our fit gives E^l'lt 7 ~ 1.45 me V. Note that 
the energy of this feature is most sensitive to J z ±. 

Moreover, we observe the following trends in the region around Yb 2 Ti 2 07: 

• as \J Z ±\ increases, all bands go up in energy (especially the two-dimensional flat one), 

• increases in J± and J±± seem to have more or less the same effect: the bands get closer, and 
this happens in particular because the energy of the top bands decreases. 

6 Mean Field Theory Calculation 

The Curie- Weiss mean field Hamiltonian obtained from Eq. (1) takes the form 



-iMBH^aT^, (20) 

i 

where W 1 is the magnetic field in the /i direction, (S'f ) is the mean field quantum thermal ex- 
pectation value of Sf, defined by (Stf) = -^Tr S? exp(— f3H MF ), where (3 is the inverse tem- 
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perature (3 = l/(/c#T), where k B is the Boltzmann constant and Z is the partition function 
Z = Tr exp(— (3H MF ). The traces are taken over the up and down spin states of every spin. 

The ground state does not enlarge the unit cell at zero temperature, and we assume this is still 
the case at non-zero temperature. Thus, we define, for every sublattice a and every axis fx the 
average magnetization 



171" 



Wt,o/> 



(21) 



which is the same for every tetrahedron t, which can be either "up" or "down". We arrive at the 
twelve consistency equations 



m„ 



7^ tanh ^^ 



2||h;f 



(22) 



where h c a ff = 2 J2b m b • Jfea — MbH • g a ; the free energy per site is 



/ = - 7 E m - • J -b ■ m 6 - — J2 ln 

* a,b ^ a 



2 cosh 



(23) 



We solve Eqs. (22) numerically, choosing the solution which minimizes the free energy in 



Eq. (23). With the field oriented along the [110] direction, we find the phase diagram shown in 



Supplementary Figure 2, which contains two phases. In zero field, the low temperature phase is 
a ferromagnet with net magnetization along any of the (100) axes. In low fields, this expands 
into a phase in which the average magnetization is not aligned with the field, but lies within a 
{001} plane. On increasing fields, a transition occurs to the high field state in which the net 
magnetization is aligned with the applied field, and no symmetries are spontaneously broken. This 
state is continuously connected to the high temperature paramagnetic phase. 
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Supplementary Figure 2: Field versus temperature phase diagram obtained from a mean-field 
analysis of the Hamiltonian of Eq. (1), for a field H parallel to the [110] direction and with the 
exchange constant values obtained with our fits, Eq. (3). The system displays a net magnetic mo- 
ment throughout the (H, T) plane. The blue region denotes a region where the total magnetization 
lies in the xy plane, and the green region is the paramagnetic phase; the two zones are separated 
by a continuous transition. In zero field, this transition takes place at T C MF = 3.2 K. At zero tem- 
perature, the amplitude of the transition field is H^ F = 1.1 T. The dark blue arrow shows the 
experimentally reported transition temperature of 240 mK: the actual transition occurs at a much 
lower temperature than that predicted by mean-field theory T C MF . 
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7 Perturbation Theory 

We show that the U(l) QSL described by Hermele et ail^il is stable to the addition of the terms in 
the symmetry-obtained Hamiltonian Eq. (4), provided all coupling constants are small with respect 
to the Ising exchange parameter J zz . In other words, we show that the U{\) QSL exists in a finite 
region of parameter space, specifically where J±, J z ± and J±± are small with respect to J zz . 

In that limit, one may apply perturbation theory. When J zz > 0, the ground state manifold 
of the unperturbed Hamiltonian is the extensively degenerate "two-in-two-out" manifold. When 



J z ± = J±± = 0, Eq. (4) can be mapped exactly onto the Hamiltonian of Ref.[T3j where the first 
non- vanishing and non-constant term in perturbation theory (above the "two-in-two-out" manifold) 
was shown to be third order in J±/J zz : 

tff g = -K Y, (S+SjS+SrS+S; + h.c.) , (24) 

{i,j,k,l,m,n}=0 

where K = —^ is a ring exchange interaction. iK" flips the spins on the "flippable" hexagons, 

J zz 5 

i.e. those with alternating up and down spins, and yields zero otherwise. This represents a 
ring move responsible for favoring the quantum superpositions of the U(l) QSL. As shown in 



Refs. 13 24 29 this ring Hamiltonian has as its ground state a U(l) QSL, whose low energy 
physics is described as the Coulomb phase of a U(l) gauge theory. This phase is locally stable to 
all perturbations^ in three dimensions, which is enough already to guarantee the persistence of the 
QSL state when the other exchange couplings are sufficiently small, i.e. when the induced terms 
in the effective Hamiltonian are much smaller than the ring coupling K. 

We can, however, go further and consider these effects explicitly in the perturbative limit, 
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which extends the discussion to the case when J z ±, J±± <C J zz but with no particular assumptions 
placed upon the magnitude of the induced terms in the effective Hamiltonian relative to K. When 
J z ± or J±± are non-zero, other "non-ring" effective Hamiltonians are allowed. In particular, the 
J z ± term gives rise to an effective third neighbour ferromagnetic Ising Hamiltonian: 

<«M>» 



3^ 2 ± 



where J@) = -j^. This term alone has six symmetry -related ordered ground states. Each of them 
consists of the choice of one of the six two-in-two-out tetrahedra, with the same pattern repeated 
on each "up" tetrahedron (these states are magnetic with their moment along the (100) directions) 
and contain no flippable hexagons, see Supplementary Figure 3. This implies that states which 





Supplementary Figure 3: One of the six ground states of H^ dMng = —^f^- J2(((i,j))) Sf S|. A 
blue sphere represents an "in" spin while a black sphere represents an "out" one. Because two of 
the four chain types contain non-alternating spins, there are no flippable hexagons in any of the six 
(equivalent) ground states. 



contain flippable hexagons represent an energy cost. We can therefore consider, H%® dhihyg as an 
analog of the Rokhsar-Kivelson term introduced by Hermele et a J.,l 13 l 25 lij"^^ = VNf, where N{ is 
the operator that counts the number of flippable hexagons. This term actually stabilizes the QSL 
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state.E322l In particular, the QSL phase grows in stability as V is increased from zero up to the point 
V — K, beyond which (for V > K), the system undergoes a first order transition to a degenerate 
set of classical unflippable states, which includes the six ordered (100) ferromagnetic ground states 
described above. 

In the case relevant here, when J(3)/K is sufficiently large, the system must undergo a transition 
to the unflippable (100) ground states. Since this model and the RK one of Ref. 13 agree on the 
phases both when J( 3 ) <C K and when J( 3 ) 3> K, it is natural to expect that the intervening phase 
diagram coincides in the two models as well. Therefore we expect that the QSL state is maximally 
stable when J(3)/K takes some value of 0(1). (For the values of the exchange constants given 

j2 J zz 

by our fits we find J(3)/K = z ^ j3 zz = 6.2, but we caution that this is probably outside the 
perturbative regime). We note in passing that the J( 3 ) exchange can also be expressed in terms of 
purely plaquette interactions, which might allow further analytical connection to the RK theory. 
We will not, however, pursue this further here. 

Inclusion of the other coupling J±±, higher order effects, and cross terms amongst the ex- 
change couplings does not lead to any new effects. Indeed, all the associated terms in the effective 
Hamiltonian assume a ferromagnetic or ring form, and can be subsumed in the above couplings. 
They are also higher order in J a /J zz , a — ±, z±, ±±. 
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8 Comparison between our Exchange Constants and those of 
Thompson et alW 

The correspondence between our effective spin-1/2 operators Sj in Eq. (1) and the full 7/2-angular- 



momentum operators J; used by Thompson et al. in Ref. 10 is given by projecting the full angular 
momentum into the ground state Kramer's doublet 

Pl/2 J 4 Pl/2 = ^^, (26) 

9j 
where P 1 / 2 is the projection operator to the ground state Kramer's doublet, and gj = 8/7 is the 

Lande factor. 

We use g xy /g z = 2.4 and g z = 1.79 for concreteness, but the results do not depend too much 

upon the details of this choice within the range of parameters found in the literature. With this 



choice, we find that the semi-formal relations between our parameters and those given in Ref. 10 



are 



Ji = 0.818^-9.08^-1.21^-2.34.7™ 

J 2 = -0.818 J lsiMg - 8.03 JL- 11.2^ + 6.16.7™ 

J 3 = 0.818^ + 4.88^ + 12.7^-2.34^ 
J 4 = 0.818^-0.523^-5.49^ + 5.62^. 

Forj; smg = 6.98 10- 2 meV, J lso = 1.90 HT 2 meV, J pi = -2.50 10" 2 meV and J DU = -2.33 lQ- 2 meV 



as calculated by Thompson et al. in Ref. 10 we obtain J\ = — 0.03 me V, J 2 = — 0.07 me V, 
J 3 = — 0.11 me V and J 4 = 0.05 me V, or, in the formulation of Eq. (4), J zz = 1.9 10~ 4 meV, 
J± = 3.0 10" 3 meV, J±± = 3.8 10" 2 meV and J z ± = -6.4 10" 2 meV. These values are rather 
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different from those found in our fits, Eq. (3) and Eq. (5), and indeed yield spin wave spectra in 
strong disagreement with experiment. In principle, one resolution of the difference could be that 
the exchange couplings are actually strongly temperature dependent, and distinctly different in the 



temperature range studied in Ref. 10 and at the low temperatures studied here. Such a change 
in exchange parameters could conceivable occur if the transition observed at 214-240mK had a 
substantial structural component. However, we do not have any a priori reason to suspect this. In 
any case, it would be very interesting to resolve the differences between the two exchange models. 
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